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We have continued our systematic investigations of the numerical simulations of lattice gauge theories in the 
dual formulation. These include: i) a more practical implementation of the quasi-local updating technique, ii) a 
thorough investigation of the sign problem, iii) issues related to the ergodicity of the various update algorithms, 
iv) a novel way of measuring conventional observables like plaquette in the dual formalism and v) investigations 
of thermalisation. While the dual formulation holds out a lot of promises in principle, there are still some ways 
to go before it can be made into an attractive alternative lattice formulation. 



1. INTRODUCTION 

We are reporting on our continuing investiga- 
tions of SU (2) lattice gauge theory in the dual for- 
mulation. A brief introduction to the formalism 
as well as to the techniques that have been devel- 
oped for numerical simulations was presented at 
LATTICE'99 ||l|. As was stated there, the par- 
tition function of the conventional LGT can be 
converted, upon using the techniques of character 
expansion and group integrations, into the dual 
form [PI (for another perspective on this issue see 
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This partition sum is defined over the dual lattice 
where {ja} hve over the links and {jb} over the 
diagonals to the plaquettes. The convention for 
the diagonals is that they connect the vertices of 
the odd sublattice. We have collectively desig- 
nated {ja},{jb} by {j}. Each cube of the dual 
lattice is seen to be spanned by 5 tetrahedra of 
which one is spanned entirely by {jh} while four 
are spanned by three {ja} and three {jb}- These 
features are illustrated in fig.l which exhibits a 
two-cube cell on the dual lattice. The b-links are 
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shown by dashed lines while the a-links are shown 
by solid lines. Each tetrahedron carries a weight 

factor which is the SU (2) 6-j symbol | ^ ^ ^ 

Periodic b.c for the original lattice is crucial for 
this construction. 

One immediate advantage of the dual approach 
is that variables are now discrete and for all prac- 
tical purposes can be represented by short inte- 
gers. Other advantages of this method have al- 
ready been commented upon in |jl[. 

1.1. Conventional Observables 

Usually when one goes to the dual description, 
observables in the original formulation become 
very complicated in the dual picture and vice 
versa (for a discussion of the Wilson loop in the 
dual picture see ||]). It is therefore curious to see 
that, in the present context, a whole class of con- 
ventional observables are just as easily measured 
in both the formulations. In the conventional for- 
mulation, the expectation value of the normalized 
plaquette is given by 
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where Np is the number of plaquettes in the lat- 
tice. In the dual picture it can be easily verified 
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that this expectation value is given by 
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where the average on the r.h.s is over the distribu- 
tion of j-values on the a-hnks of the dual lattice. 
This formula can be generalized to the case of any 
function of the plaquette variable and more im- 
portantly to the two-point correlation functions of 
plaquettes or of functions of plaquette variables. 
This allows for a direct comparison between the 
dual and the conventional pictures. 

2. UPDATES 

In we introduced two types of updating al- 
gorithms: local and quasi-local. It was already 
pointed out there that the local updates are not 
ergodic as they are unable to change integer spins 
to non-integer spins and vice versa. We illustrate 
the local update algorithm with the help of fig. 1 . 
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Figure 1. Illustrating the local updates. 



not change half-integers to integers and vice versa 
and the move is not ergodic by itself. 

In addition to the problem of ergodicity, there 
are two additional problems with this method. 
Every now and then one may face a situation 
where jmm = jmax in which case the move can 
not change the value of the link. Our code keeps 
track of this possibility through a counter irange 
and whenever this happens one skips the particu- 
lar site and moves on to the next one. The other 
problem has to do with the fact that 6j's occa- 
sionally vanish even when all the triangle rela- 
tions hold. Very rarely it so happens that every 
new value of the link is such that the product of 
6j's over the tetrahedra (four in the case of local 
updates of a-links and six in the case of b-links) 
vanishes and the new configuration is to be elimi- 
nated because it does not contribute to either the 
partition sum or to any of the expectation values. 
This is kept track of by the counter ispecial. 

2.1. Quasi-local Updates 

As already discussed in |l| , if we were to change 
a link from half-integer (henceforth referred to 
as fermionic) to integer (bosonic) or vice versa, 
in general the change proliferates throughout the 
lattice as at least two links have to be Z2-flipped 
in every triangle. What one needs are flippings 
involving the smallest cluster of links. Two cate- 
gories of such quasi-local updates were identified 
in ||l[]. At an even site of the lattice, the update 
algorithm called globeven involves the interior 
links of a cluster of 8 tetrahedra as shown in fig. 2. 
Likewise, at odd sites the corresponding update is 



The b-link shown as a dotted line is contained 
in 6 tetrahedra (two of them containing only b- 
links); it is also contained in 6 triangles again two 
of which contain only b-links. For each triangle 
containing this link, one forms jl^ax — fi + Jl ^^'^ 
Jmin — \ ji "ill where j| 2 the remaining two 
spins in the i-th triangle containing the link to 
be updated. One then forms jmax = ™'^i{jmax\ 
and jmin — iTT'O'Xiljmin} ^^"^ updatcs the link to a 
value in the range {jmin, jmax)- Clearly this does 




Figure 2. Quasi- local updates at even sites. 
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called globodd and it involves the 6 a-links and 
12 b- links in the interior of the 8-cube cluster in- 
volving 32 tetrahedra as shown in fig. 3. In both 
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Figure 3. Quasi-local updates at odd sites. 



these schemes of quasi-local updates, one keeps all 
the links on the outer surface fixed and Z2-flips 
all the interior links. The changes should be so 
made as to satisfy all the triangle inequalities. In 
jl] , the basic idea of achieving this goal was briefiy 
mentioned. This is done using variables originally 
introduced by Bargman|^ and developed by 
This will be discussed in detail in . While this 
method based on the K(agome)-variables is the 
systematic approach to the quasi-local updates, it 
is pretty involved and time consuming. A simpli- 
fied algorithm was proposed in Q which consisted 
of simply shifting all spin values on the interior 
links by i which can be seen to preserve all the 
triangle inequalities. But this can result in sub- 
stantial changes in the Boltzman weight resulting 
in abrupt changes. In the current investigations 
we modified this by following the fixed Z2-flip by 
randomly selecting a set of configurations satisfy- 
ing the triangle inequalities and then performing 
heatbath or metropolis in the subspace of these 
configurations. 

3. THE SIGN PROBLEM 

It was already emphasised in [|l| that one of the 
main drawbacks of the dual formulation is the 
lack of positivity of the Boltzman weight which 
arises out of the 6j 's not being of a particular sign. 



There we had proposed to use the trick of gener- 
ating important samples with probabilities that 
are the modulus of the weights occurring in eqn 
(1) but by absorbing the sign of the configuration 
into the observables. It was soon pointed out to 
us by Wosiek that no matter how small the lo- 
cal probability of a wrong sign, in the thermody- 
namic limit there will be an almost equal fraction 
of configurations with the two signs. This meant 
that the average value of the sign is very close 
to zero making the procedure questionable. One 
way out is to regroup the terms in eqn(l) so that 
for each group the weight is positive. But it is not 
clear how to do this. To make some progress we 
decided to adopt another strategy based on the 
fact that sign changes in the configurations come 
about through completely local sources. To see 
this, let us look at the expression for the average 
of a generic observable O: 



where {pi = '-^} with P = J2i \Pi \ ^-^d pi are the 
weights occurring in eqn(l) and are not necessar- 
ily positive. Further, is the sign of pi. Now 
if it so happens that when thermal equilibrium 
is reached the average value of O (in the usual 
sense) is insensitive to the sign of the configura- 
tion in the sense that the average over positive 
sign configurations is very close to the average 
over negative sign configurations, the numerator 
in eqn(4) is also proportional to (e)p and we can 
simply ignore the sign problem. 

We wish to argue that this is very likely to hap- 
pen for relatively local observables and also for 
the correlation functions of relatively local ob- 
servables. The reason for this is the fact that 
the source for the change of sign of the configu- 
ration is ultralocal. Combining this with the ex- 
pectation that both signs become equally likely 
as the system size increases, one is led to surmise 
that the configurational average of reasonably lo- 
cal observables as well as the correlation functions 
of reasonably local observables are insensitive to 
the sign of the configuration. So far most of our 
runs have borne out this expectation. We found 
that whenever there was a reasonable fraction of 
configurations with each sign, the averages of in- 
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ner and outer tetrahedra were indeed insensitive 
to the sign in the hmit of thermahsation (to be 
more precise, in the hmit of nearly stationary ex- 
pectation values). We are investigating this with 
other classes of observables. We also propose to 
make a systematic study of the sign problem by 
studying the way the signs fluctuate locally as 
well by studying their correlations over the lat- 
tice. 

4. METROPOLIS AND HEATBATH 

For the local updates both the metropolis and 
heatbath methods work in fairly standard ways 
because the allowed range of values for the new 
link are very small. But for the quasi-local moves 
the situation is quite different. Even in the case 
of the quasi-local updates at even sites, we are 
dealing with the space of allowed spin values on 
six links. Truncating the maximum value of spins 
at I as we did, we have 10^ values to check for 
triangle inequalities and subsequently picking a 
new element from among the list. This is a hope- 
less task. Even using the K- variable method is 
time consuming. So what we did was to apply 
the fixed Z2-flip and select the first a-link ran- 
domly from its allowed possibilities, the second 
a-link randomly from its allowed values (which 
depends on the random choice made for the first 
a-link) and so on. For metropolis algorithm, af- 
ter selecting a sextuplet of values, the new set 
was accepted or rejected according to the usual 
metropolis rule. For the heatbath method we like- 
wise selected a certain number of sextuplets be- 
fore the Z2-flip and an equal number after the flip. 
The heatbath algorithm was then applied to this 
manageable subset. At present we do not have a 
more rigorous argument to justify our heatbath 
method. For the globodd update the problems 
are even more severe and we have again adopted 
the same pragmatic approach. 

4.1. INITIAL CONFIGURATIONS 

As initial configurations we chose from among 
three classes: i) cold start - here we set all the 
a- links and b-links to j = 1. This configuration 
clearly satisfies all the triangle inequalities, ii) ex- 
treme cold start - here all a-links and b-links were 



set equal to j = 0. It is worth noting that the ex- 
treme cold start configuration is a fixed point of 
the local update algorithms, iii) hot start - here 
a-links and b-links are randomly chosen in a way 
that all triangle inequalities are satisfied. A prac- 
tical way of reaching such configurations is to first 
start with a cold start configuration and go on 
choosing links randomly in their allowed ranges. 
After sweeping the lattice sufficiently many times 
one gets to a hot start configuration. 

5. RESULTS 

We carried out our simulations on a 8^ lattice 
at various /3-values. We simulated with the ab- 
solute values of the weights given by eqn(l) and 
monitored the fluctuations in the sign of the con- 
figurations. In a few runs we noticed that the 
sign changed very few times. In those cases, the 
average values of the inner (outer) tetrahedra for 
the positive sign and negative sign configurations 
did not approach the same values even after many 
iterations. But for these cases (e)p is not close to 
zero and one can safely use eqn (4) . On the other 
hand for a majority of runs, the sign fiuctuated 
quite a bit and even when {e)p ~ 0.5, the average 
values for the two cases were very close to each 
other bearing out the surmise made earlier. This 
was so for all the classes of update algorithms, ini- 
tial states and metropolis/heatbath. We plan to 
study this issue for a larger class of observables. 

In our simulations using only the local updates, 
we found that while the results of metropolis and 
heatbath converged for given initial states, the 
results did not converge for trials with different 
initial states. We have interpreted this as due 
to the non-ergodic nature of the local updates. 
The results for the quasi-local updates were not 
so clear cut. We also ran the code with the updat- 
ing algorithms being chosen randomly; at the first 
instance we have run the code such that all up- 
dating algorithms are equally likely. Here we find 
that while the runs for the cold-start and hot- 
start converged, the results for metropolis and 
heat bath are widely divergent. We have not re- 
ally understood the reason for this. 
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6. CONCLUSION 
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Clearly a lot more work needs to be done before 
establishing the dual framework as a viable al- 
ternative for numerical simulations of Gauge and 
Spin systems. The most pressing problem is to 
gain a satisfactory understanding and control of 
the sign problem. While our surmise in this direc- 
tion is well borne out by the initial runs, we are 
yet to show its correctness in a controlled man- 
ner. One interesting possibility is to study the 
model defined by the absolute value of weights in 
its own right and we plan to undertake that study 
shortly. But we still have to understand the lack 
of convergence between the different runs before 
we can take any of the conclusions seriously. To 
this end, we are going to investigate the efficacy of 
our metropolis and heat bath methods by study- 
ing them in contexts where we can get results by 
deterministic as opposed to stochastic methods. 

Another line of attack we are contemplating is 
to look for simple models that are analytically 
tractable but which nevertheless have all the es- 
sential ingredients of the dual formulation stud- 
ied here. Indeed, such models were encountered 
in the study of analytic variational methods for 
euclidean lattice gauge theories [|| and in the an- 
alytic evaluation of the partition function of the 
unit hypercube for c? = 4 SU{2) lattice gauge 
theories ||^. We are hoping that the insights so 
gained will enable us to complete the investiga- 
tions reported here. 

There are a large class of very important prob- 
lems that belong to the same calculational scheme 
as what we have attempted here. The most no- 
table are: the d = 4 version of these results 
as obtained by Halliday and Suranyi||lO|], the 
class of models considered by Ooguri |1J and 
their non-topological descendants etc. Language 
and techniques very similar to the ones used in 
the dual picture are also becoming increasingly 
popular in treatments of quantum gravity [p^ . 
Replacing the 6j's by q-deformed 6j's would im- 
mediately relate these developments to the well 
known Turaev-Viro models and their generaliza- 
tions. Clearly it is of great importance to make 
further progress in the numerical simulations of 
theories in the dual formulation. 
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